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In many applications it is advantageous to construct effective slip boundary conditions, 
which could fully characterize flow over patterned surfaces. Here we focus on laminar 
shear flows over smooth anisotropic surfaces with arbitrary scalar slip b(y), varying in only 
one direction. We derive general expressions for eigenvalues of the effective slip-length 
tensor, and show that the transverse component is equal to a half of the longitudinal one 
£>->, with twice larger local slip, 2b(y). A remarkable corollary of this relation is that the flow 

along any direction of the ID surface can be easily determined, once the longitudinal 
component of the effective slip tensor is found from the known spatially nonuniform 
scalar slip. 
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in quantifying the effects of surface chemical heterogeneities with a different local scalar 
slip on fluid motion. In this situation it is advantageous to construct the effective slip 
boundary condition, which is applied at the imaginary smooth homogeneous, but gen- 
erally anisotropi c surface, and mimics the actual one along the tr ue heterogeneously 
slipping surface (|Vinogradova fc Belvaevll201ll iKamrin et q/.ll2010l ). Such an effective 



1. Introduction 

With recent advances in microfluidics ( Stone et al\\2004 ) , renewed interest has emerged 



condition fully characterizes the flow at the real surface (on the scale larger than the 
pattern characteristic length) and can be used to solve complex hydrodynamic problems 
without tedious calculations. 

For an anisotropic texture, the effective boundary condition generally depends on the 
direction of the flow and is a tensor, b R ff = \bff} represented by a symmetric, positive 
definite 2x2 matrix ( Bazant fc Vinogradova! 2008) 
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b cff = S 9 ( v f ^ )S_ fl , (1.1) 

diagonalized by a rotation 
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Figure 1. Schematic representation of periodic textures with scalar slip boundary conditions, 
varying in direction y. The patterns of slip boundary conditions are depicted as alternating 
stripes with piecewise constant slip lengths, but our discussion is more general and applies to 
any ID distribution of a local slip (e.g., sinusoidal, trapezoidal, and more). 



Therefore, Eq. (|l.ll) allows us to calculate an effective slip in any direction given by an 
angle 9. In other words, the general problem reduces to computing the two eigenvalues, b cS 
(9 = 0) and b^ s (9 = 7r/2), which attain the maximal and minimal directional slip lengths, 
respectively. This tensorial slip approach, based on a consideration of a 'macroscale' 
fluid motion instead of solving hydrodynamic equations at the scale of the individua l 
pattern, was supported by statistical diffusion arguments (jBazant fe Vinogradovall2008l ). 
and was recently justified fo r the case of Stokes flow over a broad class of periodic 
surfaces (JKamrin et all 120101) . 



Th e concept of an effective te nsorial slip has recently been proven by simulations (jPrieziev 



201 It ISchmieschek et aLll2012l) . and was already used to obtain simple solutions of sev- 
eral c omplex problems. It may be useful in many situati ons, such as drainage of thin 
films ( Belvaev fe V inogradova 2010fe[_ Asmolov et aLll2011 ). mixing in superhydr ophobic 
channels ( Vin ogradova & Bclyacv 20111), and electrokinetics of patterned surfaces ( Bahga et al. 
2010t iBelvaev fc Vinogradoval[2(Fllt ). However, to the best of our knowledge all these an- 
alytical solutions were obtained for a flow on alternating slip and no-slip stripes, and the 
quantitative understanding of effective slippage past other types of anisotropic surfaces 
is still challenging. 

In this paper, we study the Stokes flow past flat surfaces, where the local (scalar) slip 
length b varies only in one direction. Our focus is on the limit of a thick (compared to 
texture period) channel or a single interface, so that effective slip is a characteristics 
of a heterogeneous interface solely and does not depend on the channel thickness. We 
derive a simple universal relationship between eigenvalues of the slip-length tensor. This 
allows one to avoid tedious calculations of flows in a transverse configuration, by reducing 
the problem to an analysis of longitudinal flows, which is much easier to evaluate. Our 
results open a possibility to solve a broad class of hydrodynamic problems for ID textured 
surfaces. 



2. Theory 



2.1. General 



iideration 



We consider a creeping flow along a plane anisotropic wall, and a Cartesian coordinate 
system (x, y, z) (Fig. Q]). The origin of coordinates is placed at the fiat interface, charac- 
terized by a slip length b(y) , spatially varying in one direction, and the texture varies over 
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a period L. Our analysis is based on the limit of a thick channel or a single interface, so 
that the velocity profile sufficiently far above the surface, at a height may be considered 
as a linear shear flow. Note that our results do not apply to a thin or an arbitrary channel 
situation, where the effe ctive slip scales with the channel width ( Feuillebois et a?j|2009t 
Schmieschek et aZ.ll2012h . 



Dimensionless variables are defined by using L as a reference length scale, the shear 
rate sufficiently far above the surface, G, and the fluid kinematic viscosity, v. We seek 
the solution for the velocity profile in the form 

v = U + Ui, 

where U = ze;, I = x, y, is the undisturbed linear shear flow, e; are the unit vectors. 
The perturbation of the flow, Ui = (u,v,w) , which is caused by the presence of the 
texture and decays far from the surface at small Reynolds number Re — GL 2 jv satisfies 
dimensionless Stokes equations, 

V -111=0, (2.1) 

Vp - Am = 0, 

where p is pressure. The boundary conditions at the wall and at infinity are defined in 
the usual way 

z = 0: u lr - J 8(y)^l= J 8(v)e, J (2.2) 

oz 

w = 0, (2.3) 

z -> oo : -^- = 0. (2.4) 

oz 

where Ui r = (u, v, 0) is the velocity along the wall and /3 = b/L is the normalized slip 

length. 

A local slip length can be expanded in a Fourier series 



P (y) = X] b * ( n ) exp ( [k nv) 



n— — oo 



k„ = 2irn. 
Similarly, the solution to (|2.f|) for Ui and p has the form 

oo oo 

ui = y^ u* (n, z) exp (ikny) > P — X! P* ( n , z ) ex P(^ n y) , (2.5) 



n— — oo n— — oo 



and the Stokes equations can then be rewritten as 

V* • u* = 0, (2.6) 

VV - A*u* = 0, (2.7) 

d\ . d 2 



v*^o, ifc „,-j, A* = --e 

We have therefore reduced the problem to the system of ordinary differential equations 
(ODE), which can be now solved analytic ally. Similar strategy has been used before to 
address different hydrodynamic problems ( Asmolovll2008t iKamrin et al\\201(jt) . 



The zero mode solution represents a constant, u* (0, z) = (c x (0) , c v (0) , 0) . The eigen- 
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values of the effective slip-length tensor can be obtained as the components of u* (0, z) 



J cff 



Lc x (0) : 



Lc y (0) 



(2.8) 



Below we analyze in more details two configurations of the shear flow, where the slip 
regions are distributed parallel (6 = 0) and transverse (8 = 7r/2) to the shear flow 
direction. 

2.2. Longitudinal configuration 

For longitudinal patterns, U = ze x , so that the perturbation of the velocity has the only 
one component, Ui = (u, 0, 0) , and the system (|2.6l) - (j2.7p reduces to a single ODE: 



A*u* =0. 
Its decaying at infinity solution for non-zero modes has the form 

u* = c x (n)exp(- \k n \ z) . 



(2.9) 



(2.10) 



The boundary condition at the wall, (I2.2[) . then determines constants c x (n). Indeed, from 
(j2"3j) and (J2~TU)) we get 



z = 0: U= ^ c x (n) exp {\k n y) , — = - ^ \k n \ c x (n) exp (ifc„y) , 

n— — oo n— — oo 

so that (J2.2I) takes the form 

oo 

c z ( n ) + X! I fc m| c K (m)6* (n - m) = 6* (n) , ( 2 -ll) 

m— — oo 

where &* (n) is the Fourier coefficient of the slip length. Thus we reduce the longitudinal 
problem to the infinite linear system for c x (n) , which can be found by truncating the 
system and by using standard routines for linear systems. 

2.3. Transverse configuration 

For transverse patterns, U = ze y and Ui = (0,v,w) , the solution is usually constructed 
using the vorticity w = V x ui = (0, ,oj z ), which significantly complicates the analy- 
sis as compared to a longitudin al case (jCottin-Bizonne et al\\2004 . iPrieziev et aZ.ll2005 ; 



Belvaev fc Vinogradova! l201Qq( ). However, simp le analytical solutions can be obta ined 



directly for the Fourier coefficients of velocities ( Asmolovll2008t iKamrin et «Lll2010h . In- 
deed, equations (|2.6j) - (|2.7p for transverse stripes can be written as 

., * dw* 

ik«w H ; — 

dz 

ik n p* — A*v* 
dp* 



0, 
= 0. 



(2.12) 



dz 



- A*w* = 0. 



By excluding p* and v*, one can transforms them to a single ODE for the z-component: 

A* 2 W * = 0, (2.13) 

z = 0; z^ +oo : w* = 0. (2.14) 

A general solution of the fourth-order ODE (|2.13j) decaying at infinity can be written as 

w* = (c z +d z z)exp(- \k n \z). 
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Condition (|2.14j) then determines the coefficient, c z = 0. The velocity Fourier coefficient 
in y-direction can also be obtained from the first Eq. (|2.12j) : 

v* = c y exp (- \k n \ z) (1 - \k n \ z) , (2-15) 

with Cy = \d z /kn, so that 

w* = —ik n c y zexp(— \k n \ z) . (2-16) 

Thus the solution for a given Fourier mode involves only one unknown constant c y (n) , 
which can be found by applying boundary conditions (12.2[) . Using (|2.5[) and (12.151) . we 
then get 

OO £. oo 

z = 0: v= ^ c y (n) exp (ik n y) , — = -2 ^ \k n \ c y (n) exp (ik„y) , 

n— — oo n— — oo 

so that the slip boundary condition (J2.2D can be rewritten as 

oo 

c y (n) + 2 Y^ \k m \cy{m)b*(n-m) = b*(n). (2.17) 

m— — oo 

The system (|2.17[) is very similar to that for the longitudinal patterns, (|2.1ip . and differs 
only by the prefactor of 2. This means that the solution for c y (n) can be expressed in 
terms of coefficients c X 2 (n) for the longitudinal flow with twice larger local slip, 62 (y) = 
2b (y). Since b^ (n) = 2b* (n), we obtain from ([2~TT|) : 

OO 

c x2 {n) + 2 £ |fe m | c x2 (m) 6* (n-m) = 26* (n) . (2.18) 

m— — oc 

The left-hand sides of the linear systems (j2.17[) and f|2.18[) are identical, but the right- 
hand sides differ by the factor of 2. In what follows, 

/ s. c x2 (n) 
Cy( n ) = — 2 — ' ( ' 

Whence by using (|2.8|) we get 

bi , lHv)/L]= ^mm. (2 , 0) 

Thus, the longitudinal and transverse effective slip lengths are affinc, being related by a 
simple formula. 



3. Discussion 

In this section we compare Eq. (|2.20p with the results obtained earlier for some par- 
ticular periodic textures spatially varying in one direction y along the surface. We also 
make use of Eq. (|2.19p to prove a similarity of velocity profiles in eigendirections, again 
by giving some supporting examples from prior work. Finally, we show that our results 
are valid for arbitrary, not necessarily periodic, ID textures (varying on a characteristic 
scale L). 

3.1. Effective slip length 

During the last few decades several theoretical papers have been concerned with the flow 
past alternating (parallel or transverse) stripes characterized by piecewise constant, slip 
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lengths, b + and b~ with area fractions <p + and cj)~ = 1 — (f> + , correspondingly. Without 
loss of generality we consider below that < b~ < b + < oo. A large fraction of these 
papers deals with an ideal case of stripes with b + 3> L (perfect slip) and b~ = (no-slip). 
The formula describing th e effective slip in longitudinal direction for such a texture was 
proposed bv lPhilipl (|l972J ): 



ideal 



7T 



7T0" 1 



(3.1) 



Later lLauga fc Stond ( 2003 ) derived an expression for the transverse configuration: 



ideal 



= #■* 

2?r 



(3.2) 



which suggested that transverse and longitudinal components of the slip-length tensor 
are related as 



"ideal 



b 



ideal 



(3.3) 



This relationship is consistent with predictions of Eq. (|2.20|) . Note that Eq. (|3.3|) should be 
valid not only for regular textures, shown in Fig.[TJ but also for periodi c textures including 
sever al stripes of differing widths, e.g., for hierarchical fractal surfaces (jCottin-Bizonne et al 
20121) . 

A few authors have discussed a situation of b~ =0 and finite b + . Numerical data 



obtained by ICottin-Bizonne et all (|2004[ ) seem to satisfy Eq. (|2.20[) . and their 6JJ ff , b^ s 



asymptotically tend to t he limi ting values, bf degi and 6^ oal , when b + /L becomes large. 



Belvaev fc Vinogradoval ( 2010al ) derived analytical expressions for this case 



L 

*1 
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sec 



(=£) 



nb~ 
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In 



-.(=£)+ tan (=£) 
sec (=£) 



2tt 



2irb+ 



hi 



sec 



(=ft)+taa(=£) 



w hich are agai n in ag reement with predictions of Eq. (I2.20I) . 

Ng fc Wand j.2009) addressed the problem of effective slip lengths for stripes with 



b + 3> L and partial b , and obtained 



°eff — "ideal 



b- ± ± b~ 

~TTi "off — "ideal + ~~ ■ 



It can be seen, that this result also satisfies Eq. (|2.20j) . 

For a weakly slipping aniso tropic texture, b(x, y) <C L, th e area-averaged isotropic 
slip length has been predicted (JBelvaev fc Vinogradovall2010al ). This means that the slip 
length tensor becomes isotropic and for all in-plane directi ons, the flow a li gns w ith the 
applied shear stress. Similar conclusion has been made by iKamrin et all ( 20101 ). They 
proposed asymptotic solutions for b c g- for a weakly slipping interface, e = max \b* (n)\ <C 
1. The two-term expansions of the effective slip lengths in e can be obtained for an 
arbitrary local slip : 



b^/L = b* (0) - £ \k n \\b*(n)\\ bi s /L = b*(0)-2 £ \k n \\b*(n)\ 2 . (3.4) 
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Since 6*. (n) = 2b* (n), Eq. (|3.4[) is also consistent with the relation (|2.20|) . 

These examples fully support our predictions, but of course Eg. (12.201) is universal and 
should hold for any ID surface. 

Moreover, Eg. (12.201) still holds when the flow is unsteady. iNg fc Wand ( 201ll ) have 



considered pressure-driven oscillatory flow in a channel with walls patterned by stripes 
of b + ^> L and b~ = 0. For a thick channel they found that both real and imaginary 
parts of the effective slip lengths roughly satisfy Eq. (|3.3l) . The reason is that Eqs. (|2.10j) . 
(|2.11[) and (|2.15H2.l7|) remain valid in this case if one replaces \k n \ by \/k 2 + iw. 

3.2. Flow field 

We remark and stress that Eg. (12.191) allows one to express the entire flow field for the 
transverse configuration of patterns in terms of the longitudinal flow field, ui (y, z) — 
u[y,z,2b(y)] . Indeed, by applying the inverse Fourier transform of (J2.15J) and (|2.16|) . 
one can easily derive by using (|2.10j) . (|2.12j) and (|2.19j) : 

V= 2{ U2 + Z -d^)' W = "2^' (3 ' 5) 

du2 fi ft^ 

p=- — . (3.6) 

Whence we conclude that at the wall 

1 dv du 2 , . 

Z = 0: V= 2 U ^ d-z^^z- (3J) 

Several important conclusions follow from Q3.5p , (13.61) for alternating perfect slip (6 + 3> 
L) and no-slip stripes (b~ = 0) . In this case ui = u, and hence, at the wall the velocity 
along this ID texture is always twice that of perpendicular t o it, u = 2v. The same 
relation betwee n velocities in eigendirections has been found bv lTeo fc Khool ( 20091 ) and 



INg et al\ ( 20101 ) for the pressure-driven flow in a wide grooved superhydrophobic channel 
Analytical solutions have been obtain ed for the longitudinal velocity it and its gradient 
du/dy at the wall for the shear flow ( Sbragaglia fc Prosperettil 120071) . Eq. (|3.6p enables 



us to obtain the pressure for the transverse flow over the perfect-slip region, \y\ < <f> + /2: 

n du sin (ny) 

z = : P=--^- = — , (3.8) 

dy ^/cos 2 (Try) - cos 2 (tt</»+/2) 

The pressure over the no-slip regions, 4> + /2 > \y\ > 1/2, where du/dy = 0, is zero. The 
distribution (|3.8j) calculated for (f> + = 0.7 is shown in Fig. [2j It grows infinitely near the 
jump in b(y) (from b + ^ L to b~ = 0) at y — ±4> + /2: 

p~±sin(7r< ? !) + /2) [7rsin(7r</> + )s] _1/2 as s = (0 + /2 - \y\) -)■ +0. (3.9) 

Note that a solution for a velocity in an arbitrary direction of the undisturbed flow, 
U = z (e x cos# + By smff) , represents a superposition of solutions in eigendirections: 

Ui = ue x cos 9 + (ve y + we z ) sin 9. 

Therefore, the knowledge about the longitudinal velocity at the wall allows us to calculate 
the flow field in any direction given by an angle 9. 

3.3. Surfaces with an arbitrary ID texture 

Up to this point, our focus has been on periodic ID surface. Now, we can remove all 
periodicity requirements we used to derive Eq. (|2.19|) by applying the Fourier analysis. 
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Figure 2. Pressure distribution in the transverse flow over alternating perfect-slip, \y\ < + /2, 
and no-slip, + /2 > \y\ > 1/2, stripes. Solid and dashed lines are Eqs. (|3.8[) and (|3.9[) . respec- 
tively. 



Indeed, if ui is the solution of the longitudinal problem with a non-periodic slip length 
2Lj3(y), satisfying 

Au 2 = 0, 







"2 



oz 



2/9 (y), 



then one can verify directly that transverse solutions, Eqs. ([33]), (I3.6p . taking into account 
(|3.7|) . satisfy the Stokes equation (|2.1|) and boundary conditions (|2.2[) . ()2.3|) . Thus, our 
results are valid for any arbitrary patterned, not necessarily periodic, ID surfaces. 
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